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We use two-dimensional QCD as a toy laboratory to study 
off-forward parton distributions (OFPDs) in a covariant field 
theory. Exact expressions (to leading order in 1/Nc) are pre- 
sented for OFPDs in this model and are evaluated for some 
specific numerical examples. Special emphasis is put on com- 
paring the x > ( and x < £ regimes as well as on analyzing 
the implications for the light-cone description of form factors. 



I. INTRODUCTION 

Deeply virtual Compton scattering (DVCS) provides 
a novel tool to explore hadron structure. In contrast 
to deep-inelastic scattering (DIS), where one measures 
the imaginary part of the forward Compton amplitude 
only, DVCS allows measuring the off-forward Compton 
amplitude. From the parton point of view this implies 
that DVCS allows measuring off-forward matrix elements 
of parton correlation functions, i.e. on can access light- 
cone correlation functions of the form 

fr(x,t)= /^W(0) 7 + ^-)b)e^ + *~, (1.1) 



where — 



p 3 refers to the 



a;- 1 - = x" ± x and p + = p { 
usual light-cone components and t = q 2 = (p — p') 2 is the 
invariant momentum transfer. The "off-forwardness" (or 

skewedness) in Eq. (LI) is defined to be C = ^r-Q 



The main difference between these off-forward parton 
distributions (OFPDs) and the usual (forward ) p arton 
distribution functions is that p' =/= p in Eq. (1.1), i.e. 
since the initial state is different from the final state in 
the above matrix element, OFPDs cannot be interpreted 
as parton densities. Although it is not yet entirely clear 
what intuitive physical information is contained in the 
off- forward distributions, it is already clear that they play 
a 'dual' role between form-factors and par ton distribution 
functions : when one integrates Eq. (1.1) over x one 
obtains the matrix element of a local operator, i.e. form 
factors. On the other hand, in the forward limit ( = 
and t — 0, one recovers the usual parton densities, and 
the case ( — 0, t ^ is related to impact parameter de- 
pendent parton distributions. In a certain sense one can 



1 Note that we have chosen to follow here the conventions of 
Ref. M, which differ slightly from those in Ref. M. 



therefore interpret OFPDs as a 'parton decomposition of 
form factors', because f^(x,t) gives the contribution of 
partons with momentum fraction x to the form factor 
F(t). 

Therefore, OFPDs might be of tremendous help in 
studying the light-cone wave functions of nucleons. For 
example, in the light-cone formalism, form factors for 
q 2 < can be evaluated in the Drell-Yan-West frame 
where q + = 0. In this frame, the current is diagonal in 
Fock space and form factors can be expressed in terms 
of overlap integrals between light-cone wave functions of 
the target (see e.g. Ref. pi). Compared to the form 
factors, OFPDs give a much more detailed information 
about light-cone wave functions since they involve one 
less integration. 

One of the great difficulties in applying light-cone 
wave function based phenomenology to form factors with 
q 2 > is that one can no longer go to a frame where 
q + = and thus the current operator is not necessar- 
ily diagonal in Fock space. As a result, form factors for 
q 2 > can no longer be expressed as overlap integrals 
involving only the target's light-cone wave function, but 
one must also include contributions where the current 
first forms a meson which subsequently merges with the 
target. OFPDs provide the unique opportunity to sepa- 
rate these contributions, since the region x > £ involves 
only contributions where the current is diagonal in Fock 
space, while x < C corresponds to the regime where the 
current creates a qq pair. 

Despite the similarities between ordinary parton distri- 
butions and OFPDs in certain limiting cases, the phys- 
ical interpretation of OFPDs in general is still less ob- 
vious. This fact makes it particularly difficult to model 
OFPDs since it is not clear what properties one should 
expect and therefore toy model studies may be very use- 
ful during such a development stage. In such a situation, 
where one is studying a new observable with little known 
properties, it is often very useful to consider solvable toy 
models. 

For this purpose, we will use the 't Hooft model j|] 
(1+1 dimensional QCD for Nc — > oo) as an exactly solv- 
able toy model to study OFPDs in a relativistic field 
theory. Despite being just a toy model, QCDi + i has a 
number of interesting features, such as confinement, in 
common with QCD 3+ i. Of course, although toy models 
cannot predict what really happens in 3+1 dimensional 
QCD, they can provide useful insight on typical features 
that one should expect in a relativistic field theory: un- 
like most phenomenological models, the 't Hooft model 
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describes a covariant field theory and its solutions still re- 
flect the implicit constraints inherent to relativistic field 
theories, such as for example boost invariance and ana- 
lyticity. 

The paper is organized as follows. In Section II, form 
factors for the 't Hooft model are expressed in terms of 
light-cone wave functions, because of the close connec- 
tion between form factors and OFPDs, these form the 
basis for deriving explicit expressions for OFPDs in terms 
of light-cone wave functions in Section III. Explicit nu- 
merical examples are presented in Section IV, and the 
results are summarized in Section V. 



II. PARTON DISTRIBUTION FUNCTIONS AND 
FORM FACTORS IN TWO-DIMENSIONAL QCD 

In the 't Hooft model ||, i.e. two dimensional QCD 
in the large Nc limit, mesons are described as (non- 
interacting) bound states of a quark-antiquark pair (ab) 



and 



M 2 ^(x) 




*(») 



- vY 



(2.1) 



Here x denotes the momentum fraction carried by the 
quark in the infinite momentum (light-front) frame, g 2 = 

9 is kept fixed in the Nc — > oo limit. The eigenval- 



ues M 2 of Eq. (2.1) form a discrete spectrum and cor- 



respond to the invariant masses of meson resonances in 
this model. 

Since the 't Hooft model is formulated in a light-front 
framework, parton distribution, which have the physical 
meaning of momentum densities in the infinite momen- 
tum frame, can be easily calculated from the eigenfunc- 
tions to the 't Hooft equation. One thus finds for the 
distribution function of the quark a in state with light- 
cone wave function ^(x) M 



f a (x) = |*(s)|' 



(2.2) 



In Ref. ||, one can also find exact expressions for the 
vector form factor. For details, the reader is referred 
to this very useful paper. For example, for the matrix 
element of the + component of the vector current, which 
couples only to the quark, between two meson states one 
finds 



(m,p'b + (0)|n,p)=2P+F+„(C) 



(2.3) 



G(u;t) = / dv G(u,v;t) 



G(u,v;t)=J2 



<t>k(u)(t)k(v) 
, t-nl+ie 



(2.5) 



is the 'Gre en's function for an ad pair, i.e. the ji 2 ^ and 
4> n in Eq. (^^) are eigenvalues and eigenfunctions to the 
't Hooft equation with mi = m a . 



Here q = p — p' ', £ 



and t = q 2 . In 1+1 dimen- 



sions, there is no transverse momentum and therefore 
the invariant momentum transfer t is related to the 'off- 
forwardness' C, through energy conservation, i.e. 



Mi 



t Ml 

C i-C 



(2.6) 



In Section [II , we will provide a interpretation of OF- 
PDs in terms of light-cone time ordered diagrams. Since 
the Fock space interpretation of the various terms in form 
factors is very similar to the the terms contributing to 
OFPDs, the reader is referred to Fig. [I] for t his purpose. 

The first term on the r.h.s. in Eq. (|2.4| ) arises from 
pieces in the current operator which are diagonal in Fock 
space. Those terms can be directly evaluated from the 
overlap integral between the light-cone momentum space 
wave functions in the two particle sector. The second 
term on the r.h.s. in Eq. (2.4) corresponds to pair cre- 
ation terms in the current operator, i.e. physically it 
corresponds to the emission of a virtual meson which is 
subsequently absorbed by the external current. Naively, 
one may not have expected that such a contribution sur- 
vives the Nc — > oo limit where meson-meson vertices are 
suppressed by one power of ^== . However, the vacuum 

to meson matrix element of the current operator scales 
like y/Nc, which compensates for the suppression 
of the triple meson vertex. 

The Green's func tion G(u,v;t), in the second term on 
the r.h.s. of Eq. (2.4), describes the non-perturbative 
interaction of the qq pair emanating from the external 
current before it interacts with the target meson. Fol- 
lowing Ref. jfj), we describe this interaction by inserting 
a complete set of meson states. It should be emphasized 
that Eq. ( |2.4D is exact to leading order in i.e. it 
should have all properties that one would expect from a 
form factor in a relativistic field theo ry. For further de- 
tails, as well as a derivation of Eq. ( |2.4| ), the reader is 
referred to Ref. [01. 



where 



(2.4) 



l-Cj (y-tw) 2 



in. 



OFF FORWARD DISTRIBUTION 
FUNCTIONS 



Since Einhorn's result for the form factors (2.4) is al- 
ready expressed in terms of light-cone variables, it is a 
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straightforward exercise to isolate the off forward distri- 
bution functions. For x > £, only the piece diagonal in 
Fock space contributes and one finds (Fig. nk) 



(x > 0- (3.1) 



Note that En. ([O]) saturates the positivity constraint 
for OFPDs 0. Eq. (3T) reduces to the usual parton 
distribution [Eq. (O)] for ( = 0, i.e. 



/ c =o(a0 = /(*). (3.2) 

For a; < C, the bifocal current in Eq. ( |l . 1| ) creates a qq 
pair with momentum fraction £, i.e. only the off diagonal 
piece contributes and one finds (Fig. 0b,c) 

fax) = f b ( (x) + f c c (x) (x<Q, (3.3) 



where 



fl(x)=rjdwjdy 



and 

f c c (x) = -fjdwjdy <T 



(x < 0, 

y-( \ ^(w) r fw x 
i-t){y-wy °U'C' 



(3.4) 



y-( \ ®(y) r (w x 
i- CJ(y-w) 2 U'C 



(3.5) 



order in 

JVC 



Like Eq. (2.4), Eqs. (3.1) and (3.3) are exact to leading 



a.) 



b.) 
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FIG. 1. Light-cone time-ordered diagrams, illustrating the 
physics of the various contributions to f((x). a.) Diagrams di- 
agonal in Fock space contributing for x > b.),c.) diagrams 
off-diagonal in Fock space, which are relevant for x < (. All 
momenta are in units of p + . The larger gray blobs represent 
the (light-cone) wave function of the incoming and outgoing 
meson respectively. The two smaller blobs connected by a 
meson propagator in b.) and c.) represent the insertion of a 
complete set of meson states. 



The physical interpretation of Eq 
for x < Cj the current in Eq. (1.1) 



3.3j ) is as follows: 
creates a quark 
(with momentum xp + and an antiquark with momentum 
(C — x)p + . This qq pair can of course interact with itself 
before merging with the external meson states. Simi- 
lar to the expression for the form factor, this interaction 
among the qq pair can be described by the Green's func- 
tion G(u,v;t), which suggests the same physical picture 
of a qq pair, which forms a meson before being absorbed 
by the target meson. 



IV. NUMERICAL RESULTS 

In the numerical calculations, we used basis functions 
of the form 



Xn(x) =x r \l-xyP n {x), 



(4-1) 



where (3 is determined from the boundary condition 
7r/3 cot 7r/3 = 1 — 2pj- and P n (x) is a polynomial of n — th 
order. We found that n < 10 basis function are suffi- 
cient to achieve numerical convergence (first five figures 
of ground state masses stable and OFPDs did not show 
any visible dependence o n th e size of the basis). 

The 't Hooft equation (2T) contains two dimensionful 
parameters: the gauge coupling g and the quark mass 
to. Of course, since the 't Hooft model is only a toy 
model, there is no fundamental reason to associate any 
parameters in the model with experimental numbers in 
the real world. However, since we are trying to gain 
an intuitive understanding in general for regimes where 
quark masses are very small to 'medium sized' (u,d and 
s quarks), we adopt the following procedure to fix the 
parameters: First the coupling g is adjusted to yield the 
physical string tension in QCD, i.e. we demand 



^ 2 = 0.18(GeV)^ 



(4.2) 



yielding g = 3iOMeV. 

The quark masses are determined by demanding that 
the masses of ground state qq and qs mesons coincide 
with the masses of 7r and K mesons respectively, yielding 



Mhg 
Lip. 



(4.3) 



Numerical results for the OFPD f((x) for qq and ss 
ground state mesons are shown in Figs. @ and H respec- 
tively. Results for the first excited meson state are shown 
m Fig. | For the 7r meson in the 't Hooft model the dis- 
tribution amplitude is nearly flat. As a result, the parton 
distribution function (( = in Fig. ||) , which is just the 
square of the wave function in this model, is also nearly 
flat. For nonzero Fig. @ exhibits several interesting 
features. First of all, / c (l) = / C (C) = / c (0) = 0. Sec- 
ondly, f((x) is nearly flat for < ( < x as well as for 
C < x < 1. 
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FIG. 2. Off-forward parton distributions in a qq meson 
with light quarks (m n = UOMeV) for ( = 0, 0.3 and 0.6. 
The C = result is the conventional parton distribution. 

In the 't Hooft model, the (target-) meson wave func- 
tion is nearly flat (for small m q ) and vanishes (sharply) at 
the boundary. Since f((x) is just the overlap of the wave 
function with itself, these features are reflected in Jq(x) 
for £ < x < 1, i.e. rapid rise from zero near x = nearly 
flat for £ < x < 1 and a rapid drop near x = 1. In 3+1 
dimensional QCD one would also expect that = 

because parton distributions vanish at that point as well. 
However, since parton distributions in QCD^+i do not 



FIG. 3. Same as Fig. 
meson state. 



but for the first orbitally excited 



vanish at x = Q, one would not expect /<;(£) to vanish 
either. The physics of the regime < x < £ is entirely 
different. Here the x dependence comes from the x de- 
pendence of the distribution amplitude of the meson that 
couples to the external current. In general, that contri- 
bution is complicated since there is an infinite sum of 
mesons contributing. However, when t is near a pole in 
this sum, that particular meson will of course dominate 
and as a result the x dependence will be proportional 
to the distribution amplitude of that particular meson. 
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This is the case for small quark masses in QCTq+i! In 
two dimensions, the '7r'-meson couples not only to pseu- 
doscalar but also to vector currents and thus the V does 
contribute to the sum over meson p oles in the Green's 
function relevant for vector currents ( |2.5| ). 

In two dimensions, t and £ are not independent vari- 
ables, but are rela ted to each other (and the target mass) 
through Eq. ( |2.6| ). Since the pion mass goes to zero for 
small m q and since t is proportional to the target mass 2 
(here also the ir) we are in a situation where the ir pole 
contribution dominates in the sum over meson poles. As 
a result, the x dependence for < x < C is, to a good 
approximation, proportional to \l/ 7r (^). The coefficient 
of proportionality depends on ( in a complicated way 
[described by Eq. (3.3)]. In fact it even changes sign 
around £ ss 0.5 since there are two competing terms with 
opposite sign.0 

This close connection between OFPD for < x < ( 
and distribution amplitudes as for example been pointed 
out already in Ref. Q where it was noticed that OFPD 
evolve with the Brodsky-Lepage kernel for distribution 
amplitudes when < x < C, i.e. the UV part of OFPD 
and distribution amplitudes are the same in this regime. 

In fact, one can use this connection to measure the 
light-cone distribution amplitudes of mesons contributing 
in the sums over me son poles: the 3 + 1 dimensional 
generalization of Eq. (|3.3[) reads 



/c(M) = Et 



-Ml' 



(4.4) 



where 4> n (z) is the n-th meson's distribution amplitude, 
and c n (£,i) are some coefficients characterizing the cou- 
pling of the mesons to the target. It is therefore con- 
ceivable that, by extrapolating in t to the lowest meson 
pole, to use OFPD in the regime < x < £ as a tool 
to measure meson distribution amplitudes (e.g. the p- 
meson distribution amplitude, which is otherwise hard 
to access) 

Although /f (x) in QCD\ + i is continuous and vanishes 
for x — 0, £ and 1, these properties are not very clearly 
visible in Fig. [|due to the rapid rise of the 't Hooft wave 
function near the endpoints for small m q . These features 
are better illustrated for heavier quarks, where the wave 
function rises less rapidly near the endpoints. Results for 
'strange' quarks ^ are shown in Fig. p. 

Up to trivial kinematic factors, L dxf^(x) yields the 
form factor. For a ground state meson target, even for 
small quark masses, and for typical values of £, Fig. || 
shows that most of the contribution to this integral arises 
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FIG. 4. Same as Fig. but for a quark mass on the order 
of the strange quark mass. 

from the region x > C, i.e. from the terms that are diag- 
onal in Fock space ( |3.]| ) . This is a very surprising result 
since it should cost only very little energy to create a qq 
fluctuation which could then couple to the target. ^ Nu- 



2 As is illustrated in Fig. |l| b. and c., the gluon from the 
annihilation vertex can be absorbed either by the quark or 
the antiquark! 

3 See the discussion above on how the parameters in the 't 
Hooft equation were fixed. 



'Note that the lightest meson in QCDi+i, which becomes 



massless as m„ 



0, does contribute to the sum over meson 



poles in Eq. (3.3) 
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merically, there is of c ours e a large cancellation between 
the two terms in Eq. ( |3.3| ). These two terms correspond 
to the two possibilities that the instantaneous 'gluon', 
which emerges when the antiquark from the current gets 
absorbed, couples to the quark or antiquark in the target 
meson. Because these two terms almost seem to cancel 
each other, it is suggestive to interpret this result in terms 
of a 'screening phenomenon', indicating a 'small size 'I for 
the ground state meson in the chiral limit. However, be- 
yond this very handwaving picture, we have no intuitive 
explanation for this surprising dynamical suppression of 
the off-diagonal terms in the form factor for the ground 
state meson in this model. 



V. SUMMARY 

We have presented exact expressions for off forward 
parton distributions in the 't Hooft model (QCDi +1 with 
Nc — > oo). These expressions, together with their nu- 
merical evaluation, illustrate a number of features that 
should be generic for OFPDs in any relativistic field the- 
ory (regardless of the number of dimensions). 

Comparing ligh t-front Fock s pace ove rlap integrals for 
form factors and OFPDs Q , it is evident 

that OFPDs can be interpreted as 'parton decomposi- 
tions of form factors'. In fact, the only difference between 
form factors and OFPDs is the fact that in OFPDs the 
momentum of the 'probed quark' is not integrated over 
but rather kept fixed at momentum fraction x. Of course, 
in the 't Hooft model, where the wave functions depend 
only on one parameter, and therefore the Fock-space di- 
agonal contribution to the form factor involves only one 
integration, this means that the OFPD is given by a mere 
product of the wave function with itself. This illustrates 
that, particularly for large x, where contributions from 
higher Fock components should be suppressed also in 3+1 
dimensions, OFPDs provide much more direct informa- 
tion on the light-cone wave function of the target than 
form factors. 

While the x dependence for x > £ depends only on 
properties of the target hadron's light-cone Fock space 
wave functions, this is no longer the case for x < (. 
In the latter regime, the x dependence is governed by 
the light-cone wave functions of mesons coupling to the 
external current. In fact, Eq. (3.3) reflects that for a 



fixed value of £, each meson contributes with some in- 
dependent coefficient times its light-cone wave function. 
Although it will in general be difficult to disentangle this 



sum over meson poles into contributions from individual 
mesons, the lightest state with the appropriate quantum 
numbers should give a dominant contribution for small 
t. For example, for unpolarized OFPDs, where the 7- 
matrix structure of the probe is that of a vector current, 
one would be most sensitive to contributions from p and 
m mesons and the success of the vector meson dominance 
picture for form factors suggests that there should be lit- 
tle contamination from higher states at small t. In the 
case of polarized OFPDs, where the 7- matrix structure 
is that of an axial vector current, the dominance of the 
lowest pole (the 7r) should be even more pronounced at 
small t. 

In the 't Hooft model we found that f((x) vanishes for 
x = C, and is continuous at this point. The fact that f((x) 
vanishes for x — > £ from below comes from the fact that 
meson distribution amplitudes vanish at the end-points 
in this model and, at least if distribution amplitudes are 
close to their asymptotic form in QC-D3+1, one would 
expect the same behavior there as well. The situation 
is different for x — > £ from above, since there one would 
expect significant contributions from higher Fock compo- 
nents. In fact, f^(x) might even diverge at this point. 

In the introduction, we pointed out that OFPDs allow 
to decompose form factors into their dependence on the 
light-cone momentum of the active parton. It is very re- 
markable that, despite the fact that a vector meson dom- 
inance description for the form factor works extremely 
well for small m q only a very small portion of the 
form factor for a ground state meson comes from the 
x < ( region. Instead, for typical values of most of 
the contributions to the form factor arise from x > (, 
where the current is diagonal in Fock space and the form 
factor can be expressed as an overlap between light-cone 
wave functions. It would be very interesting to study if a 
similar suppression of the x < C contributions also takes 
place in QCZ^+i- 
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5 It is not clear which definition of 'size' really sets the scale 
here, which is why we were unable to quantify this argument 
any further. 

6 For £ > x the corresponding expressions in QCDz+i can 
for example be found in Ref. ]4|. 
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